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(54) Monitoring of activity in 
remote environments 

(57) A detector is used to monitor 
the fluid coolant pool 2 of a nuclear 
reactor and operates by detecting 
short-lived fission or activation 
products by the Cerenkov radiation 
they produce. 

The detector includes a Cerenkov 
detector head 1 connected by an 
optical fibre bundle 4 extending 
through pool 2 to a photosensitive 
detector 5. The output signal from 
detector 5 is passed to a 
comparison unit 10 and processing 
unit 1 1 for examination, evaluation 
and presentation of output signal 
data for suitable use e.g. detection 
of failed fuel or changes in density 
of the pool. An optical test unit 7 
monitors the photosensitive detector 
5. 
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SPECIFICATION 

Improvements in the monitoring of activity 
in remote environments 

5 

This invention relates to the monitoring of 
activity in remote environments, and in parti- 
cular to the monitoring of Beta particle activ- 
ity. 

1 0 It has been known for some time that when 
a Beta particle — an electron — traverses a 
transparent medium, light, called Cerenkov 
radiation, after the discoverer of the phenome- 
non, is emitted provided the velocity of the 
1 5 particle is greater than the speed of light in 
the medium. Thus the amount of the Ceren- 
kov radiation is related to the energy of the 
Beta particle (electron). 

It has been proposed to detect the leakage 
20 of fission products into a reactor water-cooling 
system by detecting the high energy electrons 
from a few selected isotopes using a counter 
which relies on the Cerenkov effect. The coun- 
ter is constituted by a photomultiplier tube 
25 viewing the water. The pulses produced by 
the photomultiplier tube are proportional to 
the Cerenkov radiation emitted, and a discrim- 
inator is set to pass only pulses produced by 
electrons with energy 5MeV or greater. This 
30 effectively separates the high energy electrons 
from both low energy electrons and photomul- 
tiplier tube noise. A choice of fission products 
which emit high energy Beta Particles in- 
cludes Br 87 (half life 55.6 sec, Beta energy 
35 8.0MeV), I 36 (86.0 sec, 6.4MeV), Rb 30 (2.7 
min, 5.7MeV), Y 94 (16.5 min, 5.4MeV), and 
Kr 88 -»Rb 88 (2.77 hrs, 5.3MeV). The employ- 
ment of the latter is particularly advantageous 
for small fuel element can ruptures, since the 
40 long half-life allows ample time for the fission 
products to diffuse into the water without 
being seriously depleted by decay. 

According to the invention, a detector for 
monitoring the fluid coolant of a nuclear reac- 
45 tor operates on the detection of short-lived 
fission or activation products by the Cerenkov 
radiation they produce, and comprises a Cer- 
enkov radiation detector head disposed in the 
fluid coolant, optical coupling means for 
50 transmitting the radiation to a photosensitive 
detector outside the fluid coolant environ- 
ment, and circuitry for receiving the output 
from the photosensitive detector for signal 
evaluation at a data receiving position. 
55 The optical coupling means may be an 
optical fibre bundle which enables non-linear 
connection between Cerenkov and photosensi- 
tive detectors. The detection applications envi- 
saged may include a failed fuel monitor by 
60 detection of short-lived (eg. half-lives from 
1 / 1 0th second to 1 50 seconds) fission pro- 
ducts in the bulk coolant by a detector situ- 
ated at a position in the bulk coolant such as 
to minimise the effects of radiation, to be 
65 employed in conjunction with more localised 



detectors, which can be of the same kind or 
of another kind such as delayed neutron de- 
tectors, situated at the fluid coolant outlets 
from individual fuel element sub-assemblies, 

70 such detection may be concerned with the 
fuel of gas-, water-, or liquid metal-cooled 
nucler reactors. Other applications include 
monitoring of the oxygen level in the sodium 
pool of a fast reactor by detection of the Beta 

75 emission from the activation product N 16 , re- 
actor power measurement by monitoring the 
Beta emission arising from the Compton scat- 
tering of prompt gamma rays from the fission 
process in the fuel, decay heat measurement 

80 by monitoring the Beta emission from the 
Compton scattering of gamma ray flux from 
fission product decay, and monitoring for 
changes in the density of the sodium in the 
pool of a fast reactor by monitoring changes 

85 in the Beta emission from changes in gamma 
ray attenuation of the sodium. 

Also relevant are the monitoring of Beta 
rays from the short-lived isotopes F 20 and 
Ne 23 . The maximum energy of these isotopes 

90 (5.4MeV and 4.4MeV respectively) is much 
greater than that for Na 24 (1 .39MeV) and thus 
a detector could be used for detecting in-core 
for variation in the F 20 and Ne 23 concentration 
and hence variation in flow in the subassem- 

95 blies — eg. due to blockage build-up. 

A typical detector of the above kind, to- 
gether with its associated circuitry will now be 
described by way of example with reference to 
the sole Figure of the accompanying drawing 

1 00 which is a diagrammatic side view partly in 
section. 

The detector includes a Cerenkov detector 
head 1 disposed in the sodium pool 2 of a 
fast reactor and including elongate shield or 
105 containment 12. Beta particles radiating in the 
pool give rise to Cerenkov radiation which is 
focused by a focusing system 3 into a beam 
which is transmitted by an optical fibre bundle 

4 extending through and out of the pool 2 to 
110a photosensitive detector head 5 with interpo- 
sition of an optical coupler 6 for a test pulse 
which can be generated for testing the equip- 
ment and employed with an optical perform- 
ance test unit 7 via lines 13, 14. The output 

1 15 signal from the photo-sensitive detector head 

5 is transmitted on line 1 5 to an amplifier 8, 
from thence to a discriminator 9 which is set 
to a range of pulse heights which it is desired 
to monitor (for example, for prompt gamma 

120 rays less than 5MeV, for fission products 
5-8MeV, and for N 16 decay greater than 
8MeV), from thence to a signal comparison 
unit 10 into which is also fed on line 16 
signals from other detectors of the same type 

1 25 when multihead monitoring is to be effected. 
A signal processing unit 1 1 is employed for 
the examination, evaluation and presentation 
of signal data in the context of the particular 
application for which the detector is intended, 

1 30 some of which are set forth earlier. There is 



also an input to the signal processing unit 1 1 
on line 1 7 in the form of an optical perform- 
ance signal from the test unit 7 for test and 
possibly also calibration purposes. It is envi- 
5 saged that this would be helpful for indication 
of any deterioration of performance due to 
radiation damage. The output signal on line 
1 8 from unit 1 1 may be used to provide a 
running visual monitor indication or may have 

10 other or additional uses, for example to pro- 
vide failed fuel detection, or other control or 
data presentation apparatus. 

It is possible with combination of signals 
from other detectors of the same type to 

15 eliminate background signals by subtraction. 
One detector can be isolated from the Beta 
activity in the sodium pool and employed to 
monitor prompt gamma flux after Compton 
scattering. 

20 Advantages include simplicity of construc- 
tion and ease of replacement, elimination of 
routing any electrical cables into the pool 
environment, and very small detector resolv- 
ing times with very high count rates. Further- 

25 more the detector head can be coated with 
various shielding materials to give differential 
attentuation of Beta particles of various ener- 
gies and hence optimisation of the detection 
of a particular source. 

30 

CLAIMS 

1 . A detector for monitoring the fluid cool- 
ant of a nuclear reactor comprising a Ceren- 
kov radiation detector head disposed in the 

35 fluid coolant, optical coupling means for 
transmitting the radiation to a photosensitive 
detector outside the fluid coolant environ- 
ment, and circuitry for receiving the output 
from the photosensitive detector for signal 

40 evaluation at a data receiving position. 

2. A detector as claimed in Claim 1 , in 
which the optical coupling means comprises 
an optical fibre bundle which establishes non- 
linear connection between Cerenkov and pho- 

45 tosensitive detectors. 

3. A detector as claimed in Claim 2, in- 
cluding an optical test unit associated with the 
photosensitive detector. 

4. A detector for monitoring the fluid cool- 
50 ant of a nuclear reactor substantially as here- 
inbefore described with reference to and as 
shown in the accompanying drawings. 
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Smith-Purcell .Radiation from Metallic and 
Dielectric Photonic Crystals 

Kazuo Ohtaka 
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Smith-Purcell (SP) radiation from a charge running near a periodically grooved solid-state 
surface has a long history of investigation since the pioneering work of Smith and Purcell(l). 
It is a kind of Cerenkov radiation, in which an evanescent photon emitted by the running 
charge turns to an observable one by the umklapp scattering by the periodic array of grooves; 
an evanescent photon originally present in the phase space outside the light line gets into the 
light cone when an umklapp shift of its momentum component takes place. 

This paper treats a theoretical study of the SP spectrum from a photonic crystal, which 
is used in place of a grating as a source of the umklapp light scattering. We use as a model 
photonic crystal a periodic array of dielectric and metallic spheres, since the system of arrayed 
spheres is treated most compactly and strictly of all in solving Maxwell's equations [2][3]. The 
result of the calculation- is expected to be applicable to any photonic crystal, inducing those 
obtained by, say, lithographic technique, with some changes added, appropriate to the system 
concerned but only quantitative in nature. Our calculation is composed of two steps, the first 
of obtaining the amplitude of the evanescent light emitted by a charge running parallel to the 
surface of the photonic crystal and the the second of treating exactly the subsequent light 
scattering by it. The SP spectrum form an one-dimensional array of spheres is treated already 
by Garcia de Abajo [4J. Our concern lies in the full three-dimensional photonic crystals. In 
a previous paper [5], we demonstrated that the band structure of a photonic crystal was 
faithfully reproduced in its SP spectrum. The dependence of the intensity on the thickness 
of the slab photonic crystal was also examined; the overall intensity increases with thickness 
to a thin system, while in a system thick enough, having roughly more than several stacking 
layers, it saturates except at the positions of the resonance peaks. The heights of the peaks 
generally continue to grow up due to the enhancement of the Q values of the resonating 
photonic bands. The saturation comes from the evanescent nature of the initial light; the 
incident light sees only a few top layers, however thick the slab may be. 

In this paper we are interested in (a) the difference of the SP spectrum between dielectric 
and metallic photonic crystals, (b) the effect of the umklapp scattering in the direction normal 
to the trajectory of charge and (c) the dependence of the spectrum with the charge velocity 

In the topic (a), it is shown that metallic photonic crystals give rise to the SP spectra 
order of magnitudes stronger than those of the dielectric case. We shall analyze this feature 
by using the Drude form of the dielectric constant for the metal spheres and calculate the 
band structure taking into account the plasmon resonances [6J. It is shown that the enhanced 
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Cerenkov Radiation in Photonic 
Crystals 

Chiyan Luo, Mihai Ibanescu, Steven G. Johnson, 
J. D. Joannopoulos* 

In a conventional material, the coherent Cerenkov radiation due to a moving 
charged particle is associated with a velocity threshold, a forward-pointing 
radiation cone, and a forward direction of emission. We describe different 
behavior for the Cerenkov radiation in a photonic crystal. In particular, this 
radiation is intrinsically coupled with transition radiation and is observable 
without any threshold. Within one particle-velocity range, we found a radiation 
pattern with a backward-pointing radiation cone. In another velocity range, 
backward-propagating Cerenkov radiation can be expected. Potential applica- 
tions include velocity-sensitive particle detection and radiation generation at 
selectable frequencies. 



When a charged particle travels inside a medi- 
um, it can drive the medium to emit coherent 
electromagnetic energy called Cerenkov radia- 
tion (CR) (/). Extensively used in particle de- 
tectors and counters (2), CR in a conventional 
material possesses three key characteristics: it 
occurs only when the particle's velocity ex- 
ceeds the medium's phase velocity, the energy 
propagates only in the forward direction, and 
there is a forward-pointing conical wavefront. 
These characteristics remain qualitatively un- 
changed even in the presence of material dis- 
persion (3-6). One possible source of unusual 
CR is in a medium with simultaneously nega- 
tive permittivity and permeability, commonly 
known as a negative-index material for its re- 
versal of Snell's law of refraction (7-72), in 
which CR is predicted to flow backward; i.e., 
opposite to the particle velocity (7). Another 
possibility exists near a periodic structure, 
where simple Bragg scattering of light can give 
rise to radiation without any velocity threshold. 
This was first confirmed by Smith and Purcell 
(13) in early experiments with electrons travel- 
ing near the surface of a metallic grating. CR 
has since been studied in one-dimensionally 
periodic multilayer stacks (14, 15), and the 
Smith-Purcell effect has been extended to near 
the surface of dielectric structures (16, 17). A 
photonic crystal (18-20), where very complex 
Bragg scattering is possible, presents a rich new 
medium for unusual photon phenomena (27- 
23). We reveal a variety of CR patterns that can 
occur in a single photonic crystal under differ- 
ent particle-velocity regimes. 

CR in a photonic crystal arises from a 
coherent excitation of its eigenmodes by the 
moving charge. Its origin lies in both the 
transition radiation, which occurs when the 

Department of Physics and Center for Materials Sci- 
ence and Engineering, Massachusetts Institute of 
Technology, Cambridge. MA 02139, USA. 
*To whom correspondence should be addressed. E- 
maik joannop@mit.edu 



the position, and / is the time. Bloch waves 
arise from the resonance of a planewave 
source with the photon dispersion relation 

i» n (k) = k-v (1) 

where co is the frequency, n is the band index, 
the wave vector k can include an arbitrary 
reciprocal lattice vector G, and co n (/r) = 
(3> n (k + G). The CR consists of modes satis- 
fying Eq. 1 (24), which can be found in k 
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charge crosses a dielectric boundary or expe- 
riences an inhomogeneous dielectric environ- 
ment, and the conventional CR, in which 
coherence is preserved throughout the medi- 
um. Both effects are incorporated in our ap- 
proach. However, unlike the Smith-Purcell 
effect, in which radiation is generated via a 
periodic grating but then propagates through 
a uniform medium, this CR is generated and 
propagates within the same crystal in the 
form of Bloch waves. The properties of these 
Bloch waves can be substantially different 
from waves in a uniform medium, leading to 
effects not previously anticipated. In one 
case, we can reverse the overall cone that 
encloses all traveling electromagnetic energy. 
In another situation, we demonstrate a back- 
ward-propagating CR behavior reminiscent 
of that predicted in negative- index materials. 
These are very general results based on direct 
solutions of Maxwell's equations and should 
find applications in particle detection and 
wave production techniques. 

For simplicity, we focus on a two-dimen- 
sional photonic crystal (Fig. I). Let a charge 
q move in the (01) direction of a square 
lattice of air holes in a dielectric, in the xz 
plane, with parameters as specified in the Fig. 
I legend. Figure 1A shows the calculated 
transverse electric (the electric field in the xz 
plane, appropriate for CR) band structure of 
this photonic crystal. We take the particle's 
motion to be in the z direction and consider a 
path where the particle does not cross dielec- 
tric interfaces (Fig. IB, inset). As a reference, 
the long-wavelength phase velocity (v c ) of 
this photonic crystal is v c = 0.44c (where c is 
the speed of light). The excited radiation can 
be determined by treating the charge as a 
source with space-time dependence 5(r-v/) = 
2 k e iA * r ~ ik ' w , that is, as a superposition of 
planewaves with different wave vectors k and 
frequencies k • v, where 5 is the Dirac delta 
function, e is the base of the natural loga- 
rithm, i is the unit of imaginary numbers, r is 
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Fig. 1. Band structure and analysis of CR in a 
photonic crystal. (A) Transverse electric band 
structure of a two-dimensional square lattice of 
air columns in a dielectric (dielectric constant e 
= 12) with column radii r = 0.4a, a being the 
lattice period. The crystal structure and the 
irreducible Brillouin zone are shown as insets. T, 
X, and M are, respectively, the center, edge 
center, and corner of the first Brillouin zone. (B) 
Method of solving Eq. 1 in k space. CR occurs 
when the k T = w/v plane (dashed line) inter- 
sects a photonic-crystal dispersion surface. Blue 
arrows indicate the group velocities of CR 
modes. (C) Method of obtaining CR cone 
shapes, The group velocities for all modes ob- 
tained in (B) form a contour, a is the cone angle 
for the mode with group velocity u, and a m 
gives the angle for the overall radiation cone 
(gray dashed lines). The angular density of 
the arrows roughly reflects the CR angular 
distribution. * 
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space by intersecting the plane <o = k • v with 
the dispersion surface u> = & n (k) (Fig. IB). 
The CR behavior in real space, however, 
must be obtained from the group velocity u = 
duldk] that is, the gradient vector to a disper- 
sion surface, which can be shown to be the 
energy transport vector for each mode (25). 
The emission angle is given by the direction 
of u. Moreover, the radiation pattern can be 
deduced from u via the group velocity con- 
tour (Fig. 1C). We plot both the charge ve- 
locity v and the group velocities // of all the 
CR modes in a velocity space (26). The 
magnitude of u is proportional to the distance 
traveled by the wavefront of the associated 
photon mode, and the magnitude of v is pro- 
portional to the distance traveled by the 
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Fig. 2. Calculated CR modes with the lowest 
frequencies for the photonic crystal of Fig. 1. 
(A) The CR emission band structure (red region) 
in the first photonic band as a function of v. 
Colored vertical lines mark representative vs for 
different CR behaviors. (B) The solution k for a 
few vs indicated by the numbers in italics. (C) 
The group velocity contours for ttje represen- 
tative values of v. 



charged panicle in the same time. Thus, the 
radiation wavefront for each mode lies on a 
"group cone" (3, 6) with its apex on the 
moving particle and a half-apex angle a, the 
angle between v - u and v. A superposition 
of all group cones gives the overall CR pat- 
tern, and the maximum angle a m of all such 
of s is the half-apex angle for the overall cone. 
Propagating Bloch modes can only exist on 
the rear side of this overall cone, whereas the 
radiation fields are evanescent on the forward 
side, and across the overall cone the radiated 
field amplitude experiences a drop. In the 
special case of a uniform material, only the 
(7 = 0 modes are excited and our approach 
yields the characteristics of conventional CR. 
For an electron traveling near a grating sur- 
face, one employs the dispersion relations of 
air along with diffraction to obtain the Smith- 
Purcell radiation. 

In the present photonic crystal, we solved 
Eq. 1 using photon bands calculated by 
plane wave expansion (27). We focused on 
the solutions with the lowest frequencies; 
higher-order modes with larger G or in higher 
bands can be analyzed similarly. Figure 2 
shows the results for the CR emission fre- 
quency (o, wave vector A, and group velocity 
u. We see that for v < < v c , the radiation 
coalesces into Smith-Purcell resonances 
around <o «* G • v. For larger v, the resonances 
merge together to form emission bands out- 
side which CR is inhibited (Fig. 2A). As v 
increases, k and u within each emission band 
are strongly influenced by the photonic band 
structure. We can identify four regimes of the 
charge velocity (Fig. 2, B and C) with four 
qualitatively different CR behaviors in a pho- 
tonic crystal: 

1) v < —0.1c corresponds to the Smith- 
Purcell regime, in which the CR arises pri- 
marily through constructive interference be- 
tween consecutive unit cells and in the first 
emission band corresponds to the first res- 
onance near k z = 2tt/<z, where a is the 
lattice period. The group velocity contour 
is approximately circular, with radius v c . 
Because v < < v c , the overall wavefront 
will be roughly circular and centered on the 
slowly moving particle. 

2) The regime v « 0.15c is unique to the 
photonic crystal: The first and higher Smith- 
Purcell resonances merge into one band, in 
which the slow photon modes near the pho- 
tonic band edge are in coexistence with the 
fast modes. As can be seen from the group 
velocity contours, in going from v = 0.1c to 
v = 0. 1 5c, the modes in the forward direction 
begin to travel slower than the charged par- 
ticle (and can even have negative group ve- 
locities), eventually producing a contour at 
v = 0.15c that winds around without enclos- 
ing v. Because in this band there are some fast 
modes in the forward direction whose u z ex- 
ceeds v, the radiation pattern here should 
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possess a backward-pointing overall cone 
(a > W2), according to the analysis of Fig. 
1CT. This is therefore the re versed-cone re- 
gime. Such a reversed cone is forbidden by 
causality (7) in a uniform passive medium. 

3) In the -0.2c < v < -0.4c regime, all 
the radiation modes in the lowest emission 
band reside in die region tt/o < k z < 2ir/a 
with u z < 0, which implies that the energy 
flows opposite to v in this band The overall 
cone is now pointing forward and, as de- 
scribed below, the radiation becomes colli- 
mated in a backward direction. CR here is 
strongly similar to the predicted behavior in a 
negative- index material. However, because 
Id - u > 0 in this regime (where k' is the 
Bloch-reduced k in the first Brillouin zone), 
the photonic crystal may be regarded as an 
effective positive-index medium (22). We 
call this the backward-flux regime. 

4) In the v > v c regime, the solution to Eq. 
1 starts from zero frequency, and the group- 
velocity contour becomes an open-ended 
curve with positive u z . The charged particle 
now travels faster than all of the excited 
modes, and constructive interference is 
achieved throughout the whole photonic crys- 
tal. This behavior is identical to CR in a 
normal dispersive medium. Hence, this is the 
normal regime. 

The CR showed less interesting transition- 
al behavior for intermediate charge velocities 
between these regimes. For example, when 
~0.4c < v < ~v c , the CR showed a mixture 
of forward and backward emissions similar to 
that depicted in Fig. 1 C. 

We confirmed our analysis by performing 
finite-difference time-domain (FDTD) simu- 
lations of radiation for a moving charge in 
this photonic crystal (28). We reproduce in 
Fig. 3 A the velocity diagrams from Fig. 2C 
and present the resulting radiation-field sim- 
ulations in Fig. 3B. To quantitatively demon- 
strate the backward radiation, we also plot in 
Fig. 3C the simulated flux through a fixed 
line perpendicular to v as a function of a). At 
the low velocity, v = 0.1c, the radiation 
shows a Smith-Purcell behavior of near-iso- 
tropic wavefront. As the velocity increases to 
v = 0.15c, the overall radiation cone indeed 
reversed as predicted. Further increasing of v 
to v = 0.3c steers the radiation to the back- 
ward direction, and if v is increased to v = 
0.6c > v c the familiar CR with a sharp for- 
ward-pointing cone as in a uniform medium 
is recovered. The angular distribution of ra- 
diation is directly visible (Fig. 3B). When 
v = 0.1c or 0.15c, the radiation is distributed 
over a wide range of emission angles without 
producing a cone of intensity maxima. For 
v = 0.3c or 0.6c, however, the CR becomes 
collimated, and a definite emission angle in 
both the forward and the backward direction 
for most of the radiation energy can be ob- 
served. In particular, the crystal-induced dis- 
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persion at v = 0.6c becomes so small that a 
strong intensity peak cone is formed and 
almost overlaps with the overall cone. A clear 
distinction between the CR behaviors for v < 
v c and v > v c is that for v > v c the measured 
flux in the z direction is positive over the 
emission band, but for v < v c certain frequen- 
cy regions appear where the flux can be 
negative (Fig. 3C). Because positive and neg- 
ative flux values can occur at the same fre- 
quency [e.g., at (u><z;/(2ttc) = 0.2], this pho- 
tonic crystal does not behave as a uniform 
negative-index medium, as noted before. An- 
other distinction lies in the field pattern and is 
similar to that reported in (J) at a phonon 
resonance frequency: The near-static, nonra- 
diating field extends beyond the overall cone 
when v < v c , whereas for v > v c , the field 
outside the cone is strictly zero. These near- 
static fields create an artificial peak around 
o> = 0 in Fig. 3C for v < v c , (as in v = 0.1c) 
which we have verified to reduce to 0 for <i> 
0 with increasing computational cell sizes 
and time steps. Finally, there are high- 
frequency radiation "tails," corresponding to 
higher-order radiation, behind the charge in 
all cases of Fig. 3B. The higher order modes 
typically have smaller group velocities and 
thus form dense forward-pointing cones of 
smaller angles, as first predicted in (6) (for 
coherently driven slow-light media). Here 
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these tails can exhibit a backward radiation 
effect as well. 

The effects presented here can be extend- 
ed to three-dimensional photonic crystals 
with little change. As a complete photonic 
band gap in three dimensions is not required, 
the crystal structure can be flexibly chosen. 
For experimental studies, an appropriate 
structure that quantitatively approximates our 
calculations could be a square lattice of air 
holes perforating a finite- height silicon slab, 
interacting with fast electrons, and operated 
near the communication wavelength. The 
same physics should apply to dielectric-in- 
air-type crystals as well. A practical issue is 
whether the radiation intensity is sufficiently 
strong to be observable for small v. Our 
numerical simulation indicates that, across a 
bandwidth of around 40%, the average radi- 
ation energy is roughly in a ratio of 1 : 1.5:4: 
20 for v/c = 0.1, 0.15, 0.3, and 0.6. Further- 
more, the radiation intensity at small v may 
even be much larger than these numbers in 
narrow bandwidths around specific frequen- 
cies, because of the unusual photon density of 
states in a photonic crystal. In conventional 
CR, for v > v c an electron can emit hundreds 
of photons per centimeter of its path. Thus, 
for the velocities studied here, the emission 
rates correspond to a range starting from 
roughly 10 and ranging up to 200 photons per 





v* 2 




0.1 0.2 

(oa/2jcc 



0.1. 012 

(oa/2nc 



0.1 02 

(Soa/2jic 



0;1 0.2 

(oa/2jic 



Fig. 3. FDTD simulation results for CR in the photonic crystal of Fig. 1. Each column represents the 
results for the value of v shown on the top. (A) Overall radiation cone shapes (dashed lines) 
deduced from the group velocity contours in Fig. 2C (B) Distribution of the radiated magnetic field 
H . Blue, white, and red represent negative, zero, and positive field values, respectively. The color 
tables are chosen separately for best illustration in each case. (C) The frequency spectrum of the 
electromagnetic flux along z through a line perpendicular to v, in arbitrary units (a.u.). 



centimeter, which should be amenable to di- 
rect experimental observation. 

A number of applications also appear pos- 
sible. Particles traveling at speeds below the 
phase-velocity threshold cannot be detected by 
conventional CR counters, and currently their 
observation relies on other devices, such as 
scintillation counters, proportional counters, or 
cloud chambers. These other devices, however, 
lack the unique advantages of strong velocity 
sensitivity and good radiation directionality 
as in conventional CR (2). With a photonic 
crystal, one should be able to achieve velocity 
selectivity and distinctive radiation patterns 
without any velocity threshold. Moreover, on 
the high-energy side, CR with a sharp radia- 
tion wavefront is possible for particles trav- 
eling through an all-air path inside a photonic 
crystal, allowing complete absence of the 
impurity scattering and random ionization 
losses inherent in a dense medium. This 
should improve the performance of present 
detectors. Finally, the CR frequency is set by 
the photonic crystal and is thus selectively 
scalable beyond optical wavelengths, open- 
ing up the possibility of flexible radiation 
sources for frequencies that are otherwise 
difficult to access. 

References and Notes 

1. L D. Landau, E. M. Liftshitz, L P. Pitaevskii, Electro- 
dynamics of Continuous Media (Pergamon. New York, 
ed. 2, 1984). 

2. J. V. Jelly, Cerenkov Radiation and Its Applications 
(Pergamon, London, "1958). 

3. I. M. Frank, NucL Jnstrvm. Methods Phys. Res. Sect. A 
248. 7 (1986). 

4. G. N. Afanasiev. V. G. Kartavenko, E. N. Magar, 
Physica 8 269, 95 (1999). 

5. T. E. Stevens. J. K. Wahlstrand. J. Kuhl, R. Merlin, 
Science 291,627 (2001). 

6. I. Carusotto, M. Artoni, C. C. L Rocca. F. Bassani, Phys. 
Rev. Lett. 87, 064801 (2001). 

7. V. G. Veselago, Sov. Phys. Usp. 10, 509 (1968). 

8. J. B. Pendry, A. J. Holden. W. J. Stewart, I. Youngs. 
Phys. Rev. Lett. 76, 4773 (1996). 

9. J. B. Pendry, A. J. Holden, D. J. Robbins, W. j. Stewart. 
IEEE Trans. Microwave Theory Tech. 47, 2075 (1999). 

10. D. R. Smith, W. J. Padilla. D. C. Vier. S. C. Nemat- 
Nasser, S. Schultz, Phys. Rev. Lett. 84. 4184 (2000). 

11. R. A. Shelby, D. R. Smith, S. Schultz. Science 292, 77 
(2001). 

12. R. A. Shelby, D. R. Smith, S. C Nemat-Nasser, S. 
Schultz. Appl. Phys. Lett. 78. 489 (2001). 

13. S. J. Smith, E. M. Purcell. Phys. Rev. 92, 1069 (1953). 

14. K. F. Casey. C. Yen, Z. A. Kaprielian, Phys. Rev. 140. 
B768 (1965). 

15. B. Lastdrager, A. Tip, J. Verhoeven, Phys. Rev. E 61. 
5767 (2000). 

16. F. J. Garcia de Abajo. Phys. Rev. Lett. 82. 2776 (1999). 

17. K. Ohtaka, S. Yamaguti, Opt. Quantum Electron. 34, 
235 (2002). 

18. E. Yablonovitch, Phys. Rev. Lett. 58, 2059 (1987). 

19. S. John. Phys. Rev. Lett. 58, 2486 (1987). 

20. J. D. Joannopoulos, R. D. Meade. J. N. Winn, Photonic 
Crystals: Molding the Flow of Light (Princeton Univ. 
Press, Princeton. NJ. 1995). 

21. M. Notomi, Phys. Rev. B 62. 10696 (2000). 

22. C. Luo, S. G. Johnson. J. D. Joannopoulos. J. B. Pendry. 
Phys. Rev. B 65. 201104 (R) (2002). 

23. C. Luo, S. G. Johnson. J. D. Joannopoulos. Appl. Phys. 
Lett. 81. 2352 (2002). 

24. The CR modes are excited with a strength propor- 
tional to the density of radiation states and qv ■ 
eVnc wnere *' is the Bloch-reduced k in the first 
Brillouin zone; C = k and e* k nC is the G-Fourier 



370 



17 JANUARY 2003 VOL 299 SCIENCE www.sciencemag.org 



Reports" 



component of the Bloch electric field, as in - 2 C 

25. A. Yariv. P. Yeh. Optical Waves in Crystals: Propaga- 
tion and Control of Laser Radiation {Wiley, New York. 
1984). chap. 6. 

26. The contour formed by u is similar to the ray surface 
in (29). though it now represents group velocities of 
different frequencies. 

27. S. C. Johnson, J. D. Joannopoulos, Opt. Express 8, 173 
(2001). 



2a The simulation uses a thick perfectly matched layer 
boundary region that overlaps 10 periods of the 
photonic crystal similar to that used in (30) and can 
absorb the Bloch waves away from a band edge. The 
moving charge is implemented as a point-like current 
density that is oriented toward the direction of mo- 
tion and whose position depends on time. 

29. M. Born, L Wolf. Principles of Optics: Electromagnetic 
Theory of Propagation, Interference and Diffraction of 
Light (Cambridge Univ. Press, New York, ed. 7. 1999). 



30. M. Koshiba. Y. Tsuji, S. Sasaki. IEEE Microwave Wireless 
Components Lett. 11. 152 (2001). 

31. Supported in part by NSF's Materials Research Sci- 
ence and Engineering Center program (grant no. 
DMR-9400334) and the Department of Defense. Of- 
fice of Naval Research. Multidisciplinary University 
Research Initiative program (grant no. N00014-01-1- 
0803). 

18 October 2002; accepted 5 December 2002 



A Reversibly Switching Surface 

Joerg Lahann, 1 Samir Mitragotri, 2 Thanh-Nga Tran, 1 
Hiroki Kaido, 1 Jagannathan Sundaram, 2 Insung S. Choi, 1 * 
Saskia Hoffer, 3 Cabor A. Somorjai, 3 Robert Langer 1 ! 

We report the design of surfaces that exhibit dynamic changes in interfacial 
properties, such as wettability, in response to an electrical potential. The change 
in wetting behavior was caused by surface-confined, single-layered molecules 
undergoing conformational transitions between a hydrophilicand a moderately 
hydrophobic state. Reversible conformational transitions were confirmed at a 
molecular level with the use of sum-frequency generation spectroscopy and at 
a macroscopic level with the use of contact angle measurements. This type of 
surface design enables amplification of molecular-level conformational tran- 
sitions to macroscopic changes in surface properties without altering the chem- 
ical identity of the surface. Such reversibly switching surfaces may open pre- 
viously unknown opportunities in interfacial engineering. 



Interfacial properties, such as wetting behav- 
ior, are defined by the molecular-level struc- 
ture of the surface (1). Diverse modification 
procedures have been used to permanently 
alter wettability (2-4). Control of wettability 
has been recently demonstrated by elegant 
methods including light-induced (5-6) and 
electrochemical surface modifications 
(7-10). These systems require chemical reac- 
tions in order to control wettability. 

We demonstrate an alternative approach 
for dynamically controlling interfacial prop- 
erties that uses conformational transitions 
(switching) of surface-confined molecules. 
Polymers have been shown to undergo con- 
formational reorientations when changed 
from one solvent to another (11) ox from one 
temperature to another (12, 13) because of 
phase transitions between a well solvated and 
a poorly solvated state (14). In contrast our 
approach maintains the system's environment 
unaltered (including solvent, electrolyte con- 
tent, pH, temperature, and pressure) while 
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using an active stimulus, such as an electrical 
potential, to trigger specific conformational 
transitions (e.g., switching from an all-trans 
to a partially gauche oriented conformation; 
see Fig. 1). Amplification of conformational 
transitions to macroscopically measurable 
changes requires synergistic molecular reori- 
entations of ordered molecules. In principle, 
this is attainable with a single-molecular lay- 
er, such as a self-assembled monolayer 
(SAM) of alkanethiolates on gold (15). How- 
ever, the dense molecular packing in SAMs 
and the. strong interactions between the meth- 
ylene groups restrict dynamic molecular mo- 
tions to the outermost atoms (16, 17). All in 
situ evidence so far indicates that applied 
electrical potentials have no effect on long- 



chain alkanethiolate monolayers on gold 
within the range of chemical stability of the 
SAM (18). In other words, conventional 
SAMs are too dense to allow conformational 
transitions and consequently do not allow for 
switching. To explore SAMs as a model sys- 
tem for switching, we must establish suffi- 
cient spatial freedom for each molecule. Once 
a low-density SAM is created, preferential 
exposure of either hydrophilic or hydropho- 
bic moieties of the SAM to the surrounding 
medium could be exploited for the switching 
of macroscopic surface properties. 

(lo-Mercapto)hexadecanoic acid (MHA) 
was chosen as a model molecule because it (i) 
self-assembles on Au(l 11) into a monolayer 
and (ii) has a hydrophobic chain capped by a 
hydrophilic carboxylate group, thus poten- 
tially facilitating changes in the overall sur- 
face properties. To create a monolayer with 
sufficient spacing between the individual 
MHA molecules, we used a strategy that 
exploits synthesis and self-assembly of a 
MHA derivative with a globular end group; 
which results in a SAM that is densely 
packed with respect to the space-filling end 
groups but shows low-density packing with 
respect to the hydrophobic chains. Subse- 
quent cleavage of the space-filling end 
groups establishes a low-density SAM of 
MHA. The spatial dimensions of the precur- 
sor molecule to be used were adapted to 
match the optimum alkanethiolate density for 
conformational rearrangements. - 

The equilibrium low-energy conformational 
state of each of the sparsely packed MHA 
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Fig. 1. Idealized representation of the transition between straight (hydrophilic) and bent (hydro- 
phobic) molecular conformations (ions and solvent molecules are not shown). The precursor 
molecule MHAE, characterized by a bulky end group and a thiol head group, was synthesized from 
MHA by introducing the (2-chlorophenyl)diphenylmethyl ester group. 
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